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ABSTRACT 

We present the results of a three-dimensional, relativistic, hydrodynamic simulation of a precessing 
jet into which a compact blob of matter is injected. A comparison of synthetic radio maps computed 
from the hydrodynamic model, taking into account the appropriate light travel time delays, with those 
obtained from observations of actual superluminal sources shows that the variability of the jet emission 
is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear 
interactions between perturbations and the underlying jet and/or the external medium. These results 
question the hydrodynamic properties inferred from observed apparent motions and radio structures, 
and reveal that shock-in-jet models may be overly simplistic. 

Subject headings: galaxies: jets - hydrodynamics - radiation mechanisms: non-thermal 



1. INTRODUCTION 

Numerical hydrodynamic simulations were initially used 
to study radiosources from their largest scales (Burns, Nor- 
man, & Clarke 1991) to the coUimation and formation of 
their associated jets (Koide et al. 2002). Due to the rela- 
tivistic nature of these sources, relativistic hydrodynamic 
simulations have become necessary to study the superlu- 
minal sources present in the nuclei of active galaxies and 
in the microquasars of our galaxy. Improving on previous 
idealized analytical calculations (Marscher & Gear 1985), 
these numerical methods are capable of studying the time- 
dependent non-linear relativistic fluid dynamics, like for 
instance the formation and propagation of shocks (Martf, 
Miiller, & Ibaiiez 1994; Duncan & Hughes 1994; Koide, 
Nishikawa, & Mutel 1996; Falle & Komissarov 1996; Aloy 
et al. 1999a). The computation of the non-thermal emis- 
sion from such hydrodynamic models provides the means 
for a direct comparison between observations and theory, 
and is hence a significant step forward towards an under- 
standing of the emission and structural variability of rela- 
tivistic jets (Gomez et al. 1995, 1997; Komissarov & Falle 
1997; Mioduszewski, Hughes, & Duncan 1997; Aloy et al. 
1999b, 2000; Agudo et al. 2001). 

Motivated by recent observations (e.g., Gomez et al. 
2000; Tingay, Preston, & Jauncey 2001; Wehrle et al. 2001) 
revealing complex jet structure and dynamics suggestive of 
being the result of multi-dimensional non-linear jet insta- 
bilities, we present a three-dimensional relativistic hydro- 
dynamic and emission simulation of a precessing jet. 

2. HYDRODYNAMIC MODEL 

Our hydrodynamic code GENESIS (Aloy et al. 1999a) 
is based on high resolution shock capturing techniques, 
best suited to describe ultra-relativistic flows with strong 
discontinuities. The code employs the method of lines 
to achieve third order accuracy both in space (using a 
Parabolic Piecewise Monotonic -PPM- inter-cell recon- 



struction) and time (by means of a third order Total Vari- 
ation Diminishing -TVD- Runge-Kutta scheme). Guided 
by observational data we simulate a jet composed of an 
ultra-relativistic plasma (adiabatic index 4/3, and specific 
energy 50c^, c being the speed of light) in pressure equilib- 
rium with an ambient atmosphere in which the jet propa- 
gates. The plasma is injected with a bulk Lorentz factor 
of 6, and is 10^ times lighter than the environment. We 
enforce a twofold precession of the jet by imposing helical 
perturbations on the injection velocity with amplitudes C,i 
(C2) = 0.035 (0.005), and periodicities n (ra) = 400 ijjc 
(25 i?b/c), where Rb is the beam radius (the definitions of 
Ci_2 and Ti^2 are the same as in Aloy et al. (1999b)). The 
first, long wavelength precession (Ai = 400 i?b) causes a 
slight bending of the jet (within the computational do- 
main of longitudinal size 90i?f,), while the second preces- 
sion, with a much shorter wavelength (A2 = 25i?b) and 
smaller amplitude than the first one, gives rise to a small 
helical modulation of the beam. Both structural charac- 
teristics are observed in many astrophysical sources (e.g., 
3C 120; Gomez et al. 2000). 

According to the recent observations of Marscher et al. 
(2002) radio components are generated by the injection 
into the jet of material coming from the inner accretion 
disk. Motivated by these observations, we perturb the jet 
by injecting a blob of relativistic particles through the jet 
inlet during 0.8 Rb/c in the frame attached to the source 
(LAB-frame). The blob, initially a cylindrical region of 
thickness 0.8 Rb and width 1.0 i?b, is four times denser 
than the jet plasma, but has the same velocity and in- 
ternal energy per particle. As the injected perturbation 
travels downstream it spreads asymmetrically along the 
beam and splits into two regions (Fig. 1) that have distinct 
observational signatures (§ 3). The fact that injected jet 
perturbations experience substantial evolution is general, 
i.e., is not specific to the perturbation that we have chosen. 
Indeed, a similar hydrodynamic evolution should be 
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Fig. 1. — Three dimensional ray cast view of the simulated jet in 
the LAB-frame. The image is produced by ray tracing the Lorentz 
factor and the pressure, assigning an opacity to each volume element 
proportional to the magnitude of each variable. The computational 
domain is a cuboid of physical size 14 X 14 X 90 beam radii cov- 
ered by a uniform Cartesian {X,Y,Z) grid of 168x168x1080 zones. 
Those parts of the fluid having Lorentz factors larger than the bulk 
Lorentz factor of the beam in the injection nozzle {Wi, = 6) show 
up as pink regions. Strong yellowish green colors encode the energy 
density times the beam mass fraction {X = Pb/p, pb and P being the 
densities of beam and of the fluid, respectively) and diffuse green- 
ish shades the energy density itself. In both cases, the intensity is 
a measure of the corresponding quantity's value (brighter meaning 
larger). The two insets are zooms of the perturbation, the bottom 
one also showing regions A and B defined in the text. 

expected for any perturbation that either does not match 
the shock jump conditions at its boundaries or propagates 
along a non-uniform beam. An ideahzed ID model of the 
perturbation allows us to explain the basic features of its 
evolution. The original state disappears after « 15 Rb/c 
and displays a noticeable split after 30 Rb/c and beyond 
in the LAB-frame (Fig. 2). In the 3D simulation the evo- 
lution of the perturbation is also affected by its interac- 
tion with the external medium and changing beam con- 
ditions downstream. The front region (A) of the original 
blob shows the largest Lorentz factor (> 9) and relatively 
small energy density, while the back region (B) possesses 
the largest energy density but the smallest Lorentz factor 
(« 4.5). 

Initially the perturbation is ballistic as a consequence of 
its high relativistic density (« 9.7 times that of the exter- 
nal medium) and moves in a straight trajectory until 
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Fig. 2. — Snapshots of the ID evolution of the energy density 
(solid line) and the Lorentz factor (divided by 200; dashed line) of 
a square component in a uniformly moving background medium as 
seen in the O-frame for a viewing angle of 15° (top panels) and in 
the LAB-frame (bottom panels). Physical parameters of both back- 
ground and perturbation are the same as in the 3D simulation. The 
quasisteady state (rightmost panels) is very similar in both frames, 
but due to the time dilation effect the apparent size is much larger 
in the O-frame. 

it impacts on the bent walls of the beam, between times 
20 and AO Rb/c. This causes an increase of the internal 
energy density and a reduction of the Lorentz factor in 
the region of maximum shear (Fig. 3b, 3c). The heating 
is a consequence of the compression of the perturbation, 
which is enhanced by the presence of a standing shock in 
the beam (Fig. 3a). The interaction also leads to the for- 
mation of a conical shaped bow shock that continues along 
with the perturbation expanding at ~ 0.2 c, which is seen 
in Fig. 3a as a strong increase of the energy density emerg- 
ing from the component and forming a small angle with 
the beam. In the bottom part of the component the bow 
shock is weaker. The effect of the collision against the 
edge of the beam is relatively small because of the very 
small impact angle of ^ 5°. A 2D analytic modeling of 
the impact shows the generation of a rarefaction instead 
of a shock for impact angles smaller than ~ 35°. [This is 
a genuine relativistic effect (Pons, Martf, & Miiller 2000; 
Rezolla & Zanotti 2002).] After ^ 55 Rb/c the perturba- 
tion is no longer ballistic and no longer covers the whole 
beam's width. 



3. EMISSION 

We compute the synchrotron emission from the hydro- 
dynamic model assuming that the magnetic energy density 
is proportional to the particle energy density and that this 
proportionality remains constant throughout the whole 
computational domain (as in Gomez et al. 1995). We also 
assume that magnetic field is dynamically negligible and 
we account for the appropriate relativistic effects (includ- 
ing light travel time delays as in Gomez et al. 1997), aimed 
to obtain synthetic radio maps (Fig. 4) that can be com- 
pared with observations of superluminal sources. To relate 
emission features in the observer's frame (O-frame; 
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Fig. 3. — Longitudinal slice through part of the jet showing the 
logarithm of the energy density, e, (in arbitrary units) when the 
perturbation passes through the first recoUimation shock (a). The 
bottom panels show the logarithm of the energy density (b; with the 
same scale as the panel a) and the bulk Lorentz factor, W, (c) for 
a transverse slice at z=11.5 Ri, (vertical white line in panel a). In 
all the panels the white lines are isocontours of beam mass fraction 
(see caption of Fig. 1) 0.5 roughly separating the beam from the 
external medium, the vertical axis provides the scale in Ri,, and the 
crosses mark the position of the jet axis if it were axisymmetric. 
The spiral structure present in the ambient medium (b) is caused 
by the precession of the jet. 

Fig. 4) with their hydrodynamic counterparts, we provide 
a space-time diagram (Fig. 5). 

Figure 4 shows the appearance of a large region of in- 
creased emission, associated with the imposed perturba- 
tion, passing by several standing knots (SI, S2, S3) caused 
by jet recoUimation shocks. The large extension of this 
new emitting region is caused by the light travel time de- 
lays between the front and back of the perturbation, which 
stretch it by a factor (1 — (3cos0)~^ {f3 and 9 being the 
velocity -in units of c- of propagation and the viewing 
angle, respectively). As a consequence, the structure of 
the perturbation is magnified leading to brightness distri- 
bution variations within the component in the 0-frame. 
This may have significant implications when interpreting 
the observations of supcrluminal sources. First of all, the 
quite common practice of identifying brightness peaks in 
radio maps with components may be misleading, as the 
component selected in this way may only be one of many 
brightness features caused by a single perturbation and 
may not be related to any physical structure on its own. 
Moreover, a component may result from different regions 
of the perturbation having different hydrodynamic proper- 
ties (Fig. 4). For instance, identifying the brightness peak 
in the newly ejected emitting region in Fig. 4 as compo- 
nent M leads to variations in its apparent motion from 1.4 
to 4.2 c. The LAB-frame pattern Lorentz factors of the 
perturbation inferred from these values are 2 to 4.3, which 
cannot be associated to any flow or pattern speeds in the 
jet. On one hand, they are not related to the motion of 
any fluid element as the bulk Lorentz factors (Fig. 2) in 
the perturbation dynamically vary from the initial value 
of 6 (uniform within the perturbation) to values ranging 
from 4 (region B) to 9 (region A). On the other hand, the 




Fig. 4. — Time sequence of the simulated radio maps (total in- 
tensity in arbitrary units using a square root brightness scale) as 
seen in the O-frame. The epoch is shown at the right of each snap- 
shot. The maps are computed for a viewing angle of 15° and an 
optically thin frequency of 22 GHz. Yellow underbrackets indicate 
the extension of the imposed hydrodynamic perturbation (For the 
last two epochs arrows indicate that the front of the perturbation is 
undetermined.) The left and right bounds are associated to the red 
triangles and asterisks in the inset of Fig. 5, respectively. Knots in 
the emission pattern (SI, S2, S3, M and HI) are marked with dif- 
ferent symbols, and a white arrow gives the position of component 
HI, as determined from Fig. 5. 

pattern motions in the perturbation that correspond to the 
shocks confining regions A and B have Lorentz factors of 7 
and 4, respectively, which differ significantly from the esti- 
mated values. The proper motion of M is almost constant 
during the last four epochs and equals 1.7 c, correspond- 
ing to a LAB-frame pattern Lorentz factor of 2.2, which is 
only slightly smaller than the Lorentz factor of the back 
shock of region B during these epochs (w 2.7; Fig. 5). 

The internal brightness distribution of the supcrluminal 
component M also changes due to the interaction of the 
perturbation with the ambient medium and the jet. The 
azimuthal position of the region of maximum interaction 
(corresponding to the maximum energy density in region 
B) varies with time because the perturbation moves (dur- 
ing its ballistic epoch) along a straight trajectory while the 
beam is helical. The resulting interaction, together with 
the squeezing of the perturbation when it passes through 
standing knots in the jet, leads to transverse motions of the 
component's emission peak. A similar behavior has been 
observed in the jet of 3C 120 (Gomez et al. 2001), support- 
ing the idea of a precessing jet in this source and of inter- 
actions between components and the external medium. 
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Fig. 5. — Space-time diagram showing the evolution of the energy 
density along the jet. Every vertical pixel line is a snapshot of the 
jet's longitudinal structure averaged across azimuthal slices. The 
trajectories of two "elbows" (i.e., local bends of the helix that have 
evolved non-linearly and form kinks in the beam where the energy 
density and the density are large and leave a trace in the space-time 
diagram and in the radio maps of Fig. 4) of the secondary helix 
(helix 1 and helix 2; white dashed lines) of the helical jet propagate 
with pattern Lorentz factors between 2.0-2.6. The perturbation has 
average bulk Lorentz factors of 3.7 and 7.1 for the back (B) and 
front (A) regions, respectively. The back shock of region B decel- 
erates between 40 and 55 -Ri,/c reducing its pattern Lorentz factor 
to 2.7. The area between the two black lines (isocrones in the O- 
frame corresponding to times ORb/c and 5Rb/c) covers the whole 
time evolution of Fig. 4. The propagation of the component HI (= 
helix 1) and of the main perturbation (red asterisks: front shock of 
region A; red triangles: back shock of region B) are illustrated both 
in the LAB- and in the O-frame (inset). 

While region B generates the most significant features 
in the radio maps, no observable knot is produced by re- 
gion A because Lorentz factors are higher in this region 
and the radiation is thus beamed into a cone smaller than 
the viewing angle. A larger viewing angle would beam 
region A and the resulting component would present a 
faster superluminal motion. Hence, the jet viewing angle 
may introduce a selection effect determining which parts 
of the perturbation appear as superluminal components 



(with different apparent speeds) in observed radio maps. 

In contrast to M, the component HI suddenly appears 
at some distance further downstream in the jet (Fig. 4). 
This pop-up component (PUC) results from the non-linear 
evolution of the secondary helical motion of the beam, and 
hence does not correspond to any bulk fluid motion. Its 
apparent speed is 1.4 c, which agrees with its LAB-frame 
pattern Lorentz factor of 2.0 (Fig. 5). PUCs occur when 
a relativistic perturbation (observable or not) causes en- 
hanced beamed emission which suddenly becomes observ- 
able. In our model the enhanced emission is produced by 
shock heated matter at the edge of region A of the per- 
turbation, and the beaming is produced by the relativistic 
flow speeds in that region. The pop-up effect arises when 
the emitting region moving along the helical beam turns 
towards the observer, and the edge of the beamed radiation 
suddenly crosses the line of sight. PUCs appear to have 
been observed in actual sources (0420-014; Zhou et al. 
2000) and provide a good example of a coupling between 
phase (helical) and bulk (the fluid in region A) motions 
leading to the formation of new components. 

4. CONCLUSIONS 

Our simulation shows that the non-linear hydrodynamic 
evolution of perturbations can determine the observed ra- 
dio emission properties of superluminal sources, and that 
the interpretation of observed radio maps is error-prone 
when naively associating single shocks to superluminal 
components. Indeed, if the radio components of actual 
sources are correctly represented by those studied in our 
model, most observable features should not be related to 
fluid bulk motions, but instead to a complex combination 
of bulk and phase motions, viewing angle selection effects, 
and non-linear interactions between perturbations and the 
underlying jet and/or ambient medium. Further simula- 
tions spanning a wide range of the relevant parameters are 
underway. 
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